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Abstract 

Atomic  force  microscopes  employ  stacked  or  cylindrical  piezoceramic  actuators  to  achieve  sub-angstrom  reso¬ 
lution.  While  these  devices  produce  excellent  set-point  accuracy,  they  exhibit  hysteresis  and  constitutive  nonlin¬ 
earities  even  at  low  drive  levels.  Feedback  mechanisms  can  mitigate  the  deleterious  effects  of  these  nonlinearities 
for  low  frequency  operation  but  such  techniques  fail  at  higher  frequencies  due  to  increased  noise  to  signal  ratios. 
In  this  paper,  we  quantify  the  hysteresis  and  constitutive  nonlinearities  through  a  Preisach  model.  As  illustrated 
through  a  comparison  with  experimental  data,  this  provides  a  characterization  which  is  sufficiently  accurate  for 
inclusion  as  an  inverse  compensator  in  various  control  designs. 
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1.  Introduction 

Atomic  force  microscopes  are  a  subclass  of  scanning  probe  microscopes  which  have  redefined  diagnostic 
capabilities  in  fields  ranging  from  material  science  to  cell  biology.  To  illustrate  the  construction,  capabilities  and 
and  current  research  directions  concerning  atomic  force  microscopy,  consider  the  prototypical  design  depicted  in 
Figure  la.  Longitudinal  and  transverse  actuation  is  provided  by  either  a  piezoceramic  shell  or  a  stage  forced 
by  a  stacked  actuator  as  depicted  in  Figure  2.  As  the  sample  is  moved  laterally,  displacements  in  an  adjacent 
microcantilever  are  sensed  by  monitoring  changes  in  the  reflected  laser  path.  Corresponding  forces  are  determined 
from  Hooke’s  law  and  a  feedback  law  is  used  to  determine  the  distance  through  which  the  sample  should  be  raised 
or  lowered  to  maintain  constant  forces.  After  one  sweep,  this  produces  a  profile  of  the  sample  as  illustrated  in 
Figure  lb.  Details  regarding  the  construction  and  applications  utilizing  atomic  force  microscopes  and  scanning 
tunneling  microscopes  (STM)  can  be  found  in  [8]. 

The  accuracy  of  the  atomic  force  microscope  is  naturally  correlated  with  the  accuracy  of  the  stages  employed 
for  lateral  ( x-y )  and  transverse  or  vertical  ( z )  transduction.  Two  configurations  which  are  currently  employed 
are  stacked  actuators  utilizing  motion  and  cylindrical  actuators  poled  for  cLi  motion.  The  stacked  actuator 
offers  the  advantage  of  simplified  construction  whereas  the  cylindrical  configuration  offers  better  isolation  from 
exogenous  vibrations  in  addition  to  small  stage  construction.  While  both  configurations  provide  highly  repeatable 
and  accurate  set  point  placement,  the  relations  between  input  voltages  and  generated  strains  or  displacements 
exhibit  hysteresis  and  constitutive  nonlinearities  as  illustrated  in  Figure  3  with  data  collected  from  an  AFM 
employing  a  stacked  actuator.  At  low  frequencies,  feedback  mechanisms  can  adequately  attenuate  these  effects 
thus  leading  to  the  success  of  the  technologies.  However,  at  the  higher  frequencies  required  for  current  and  future 
applications,  two  phenomena  degrade  the  accuracy  achieved  by  present  control  techniques:  (i)  The  piezoceramic 
materials  exhibit  increased  hysteresis  and  (ii)  The  noise  to  signal  ratio  in  the  sensing  component  of  the  device 
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Figure  1.  (a)  Configuration  of  a  prototypical  AFM;  (b)  Surface  image  determined  after  one  lateral  sweep. 


increases  to  the  point  where  feedback  mechanisms  are  amplifying  noise  rather  than  attenuating  unmodeled 
dynamics. 

A  number  of  techniques  have  been  developed  to  quantify  hysteresis  in  piezoceramic  materials  including  phe¬ 
nomenological,  macroscopic  models  [6,  7,  15],  quasi-macroscopic  energy-based  models  [12,  13,  16]  and  microscopic 
models.  In  this  application,  the  requirement  of  real-time  implementation  precludes  the  use  of  microscopic  models 
due  to  the  large  number  of  required  parameters.  We  focus  here  on  the  development  of  an  appropriate  Preisach 
model  due  to  the  generality  of  the  formulation  and  the  rich  mathematical  framework  for  the  theory.  The  accu¬ 
racy  of  this  model  is  comparable  to  that  obtained  with  an  energy-based  formulation  [14]  while  providing  a  more 
general  framework.  However,  the  disadvantage  of  the  technique  lies  in  the  large  number  of  required  parameters 
and  difficulties  posed  by  operating  conditions  having  potential  temperature  variations  and  broadband  require¬ 
ments  (see  [5,  10]  for  general  discussion  regarding  techniques  for  accommodating  temperature  and  frequency 
dependence  in  Preisach  operators). 

Appropriate  constitutive  relations  are  summarized  in  Section  2  and  a  Preisach  model  quantifying  the  hys- 
teretic  relation  between  input  fields  and  the  polarization  generated  in  the  piezoceramic  positioning  elements  is 
presented.  The  performance  of  the  model  is  illustrated  in  Section  3  through  a  comparison  with  experimental 
data  plotted  in  Figure  3. 


Figure  2.  Actuator  configurations  employed  for  positioning  in  an  atomic  force  microscope,  (a)  Rectangular 
stacked  actuator;  (b)  Cylindrical  transducer. 
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Figure  3.  Relation  between  the  input  field  E  and  displacements  generated  by  a  stacked  PZT 
mechanism  in  an  AFM. 


positioning 


2.  Constitutive  Models 

We  consider  the  quantification  of  the  nonlinear  and  hysteretic  relation  between  input  fields  E  and  strains  e 
generated  in  the  piezoceramic  actuator  in  two  steps;  (i)  the  quantification  of  the  relation  between  E  and  the 
polarization  P  and  (ii)  the  characterization  of  the  relation  between  P  and  e.  Based  on  experimental  data  in  [13], 
we  make  the  assumption  that  for  the  drive  levels  under  consideration,  hysteresis  and  constitutive  nonlinearities 
are  primarily  manifested  in  the  relation  between  E  and  P.  It  is  this  component  that  we  characterize  through  a 
Preisach  model.  We  further  assume  that  the  relation  between  P  and  e  is  approximately  linear  as  substantiated 
by  classical  references  [4,  9].  Experiments  to  validate  this  latter  assumption  are  currently  being  performed. 

To  simplify  the  discussion,  we  focus  on  the  development  of  a  constitutive  model  for  the  stacked  actuator 
depicted  in  Figure  2a.  The  extension  of  the  1-D  constitutive  relations  to  the  cylindrical  geometry  is  addressed  in 
[14]  and  a  Preisach  model  for  that  actuator  configuration  can  be  developed  by  employing  the  Preisach  relations 
summarized  in  Section  2.2  in  the  2-D  constitutive  expressions  in  [14]. 

2.1.  Constitutive  Relations 

To  motivate  an  appropriate  form  for  the  1-D  constitutive  relations  for  piezoceramic  materials,  consider  first 
the  case  of  an  idealized  material  which  exhibits  negligible  hysteresis,  constitutive  nonlinearities,  or  structural 
damping.  Let  a  denote  the  stress  in  the  material  and  let  sp  ,cp  =  1  / sp  respectively  denote  the  Young’s  modulus 
and  compliance  at  constant  polarization.  As  illustrated  in  [9],  consideration  of  the  elastic  Gibbs  free  energy  for 
piezoelectric  materials  then  yields  the  linear  constitutive  relations 

e  =  spa  +  (3P 
E  =  -ba  +  xaP 


or  equivalently 

a  =  cpe  —  cp/3P 

(!) 

P  =  X£o  E  +  7  a. 

Here  6,7  and  xa  =  l/(X£o)  are  piezoelectric  and  electric  coupling  coefficients,  y  denotes  the  electric  susceptibility, 
and  £0  is  the  permittivity  of  free  space.  To  express  the  actuator  relation  in  (1)  in  terms  of  input  voltages  V,  the 
linear  polarizability  relation  P  =  x£oP  can  be  invoked  along  with  the  approximate  formulation  V  =  Eh,  where 


h  is  the  material  thickness,  to  obtain 


a  —  cpe  —  cpdV/h 
P  =  X£o  E  +  7  cr. 


The  specific  nature  of  the  piezoelectric  coefficient  d  =  Px£o  depends  on  the  configuration  of  the  actuator  and 
manner  through  which  it  is  poled.  The  stacked  actuator  utilizes  a  d33  effect  in  which  fields  in  the  3  direction 
generate  displacements  in  the  3  direction.  Alternatively,  transverse  displacements  in  the  cylindrical  transducer 
are  generated  by  the  effect  in  which  fields  through  the  shell  thickness  generate  longitudinal  (1-direction) 
displacements. 

To  incorporate  internal  damping,  saturation  nonlinearities,  and  hysteresis,  we  generalize  (1)  to  obtain 

cr  =  cpe  +  cue  —  cp/3P(E,  a) 

P  =  P(E,  a) 

where  cd  is  the  Kelvin-Voigt  damping  parameter  and  T  quantifies  the  hysteresis  and  constitutive  nonlinearities 
inherent  to  the  materials. 

One  technique  for  specifying  T  is  based  on  the  domain  wall  models  detailed  in  [12,  13]  and  an  AFM  model 
based  on  this  theory  is  provided  in  [14] .  Alternatively,  one  can  specify  T  through  a  Preisach  characterization  as 
discussed  in  the  next  subsection. 


2.2.  Preisach  Model 


To  quantify  the  nonlinear  hysteretic  relation  between  the  input  field  and  generated  polarization,  we  consider 
a  general  Preisach  expansion  in  which  the  hysteresis  is  represented  as  an  expansion  based  on  piecewise  constant 
kernels.  More  general  discussion  regarding  Preisach  models  for  piezoceramic  materials  can  be  found  in  [6,  7,  15] 
while  analogous  theory  for  magnetostrictive  and  shape  memory  compounds  can  be  found  in  [1,  2,  3,  5,  10,  17]. 

To  motivate  a  general  Preisach  kernel,  we  consider  first  a  single  relay  operator  ks.  This  kernel  is  characterized 
in  terms  of  crossing  times  r(<)  defined  by 


r(t)  =  {r)  G  (0,  t\  |  u{rf)  =  s3  or  u(r, ])  =  s2} 

where  u  =  E  denotes  the  input  function  and  s  =  (si,  s2)  are  points  in  the  Preisach  half  plane 

S  =  |s  £  IR2  |  s  =  (si,s2),si  <  s2}  . 

The  values  si,s2  are  threshold  values  for  the  multivalued  kernel  as  reflected  in  the  definition 


[ks(u,  0](i)  =  <  -1 


if 

r(t) 

-i 

if 

r(t) 

+i 

if 

T(t) 

if  r(t)  7^  0  and  u(maxr(t))  =  si 

if  r(t)  0  and  u(maxr(t))  =  s2  . 


A  depiction  of  this  kernel  is  given  in  [15].  The  starting  value 

—1  if  u(0)  <  Si 

[M«,0](°)  =  {  £  if  Si  <  «(0)  <  s2 
+1  if  u(0)  >  s2 


defines  the  initial  state  of  the  kernel  in  terms  of  the  parameter  £  £  {—1,1}. 

The  output  remains  on  a  branch  until  a  threshold  is  reached  in  the  monotonically  increasing  input  u.  At 
that  point,  the  output  jumps  to  the  other  saturation  value  and  remains  there  until  the  other  threshold  value  is 
reached.  For  example,  an  output  response  starting  with  a  value  of  —1  will  retain  that  value  until  u(t)  reaches 
s2.  The  output  then  jumps  to  +1  until  the  threshold  value  of  si  is  reached. 


The  classical  Preisach  operators  are  then  defined  in  terms  of  parallel  collections  of  these  single  relay  operators. 
To  this  end,  we  let  A4  denote  the  set  of  all  finite,  signed  Borel  measures  on  S  and  let  £  be  a  Borel  measurable 
function  mapping  S  — >  {—1, 1}.  For  u  €  C[0,T]  and  /i  G  M,  the  Preisach  operator  is  defined  by 

[P^u,  £)](*)  =  J\ks{u,t>(s))\{t)dn{s). 

The  goal  in  the  parameter  identification  problem  is  to  estimate  /z  so  that  a  model  response  “fits”  experimental 
data  in  a  least  squares  sense.  Details  regarding  the  convergence  property  of  the  operator  with  regard  to  both 
parameters  and  time,  and  extensions  to  Krasnoselskii-Pokrovskii  kernel  which  avoid  the  convergence  difficulties 
are  summarized  in  [15]. 

To  numerically  implement  the  Preisach  model,  it  is  necessary  to  approximate  both  the  operator  and  measure. 
As  detailed  in  [1,  15],  the  discrete  operator  is  given  by 

m 

[Pm(u,Q](t)  =  ^2 l^rii  (u,  )]  (t) ani  (3) 

i= 1 

where  n.j  are  nodes  in  the  restricted  Preisach  plane  S a  =  {(si,  S'2)  | <  Si  <  S2  <  smax}  and  aHi  are  weights 
in  the  expansion 

m 

l+rn  —  ^  (4) 

i=l 

used  to  approximate  the  measure  /z  by  a  linear  combination  of  Dirac  measures  with  atoms  at  n,;.  Details 
regarding  the  convergence  properties  of  the  approximate  operator  and  measures  are  provided  in  [1]  while  an 
inversion  technique  utilizing  a  discrete  Preisach  operator  of  this  form  is  experimentally  implemented  in  [17].  The 
performance  of  the  operator  (3)  for  quantifying  the  hysteresis  and  constitutive  nonlinearities  in  the  piezoceramic 
actuators  employed  in  the  AFM  is  illustrated  in  the  next  section. 


3.  Model  Validation 

To  illustrate  the  predictive  capabilities  of  the  model,  we  consider  the  characterization  of  hysteretic  data 
collected  from  the  stacked  actuator  depicted  in  Figure  2a.  This  data  was  collected  at  low  frequency  and  hence 
represents  a  quasistatic  application  of  both  the  device  and  model.  As  illustrated  in  Figure  4,  the  data  is  both 
hysteretic  and  exhibits  certain  saturation  nonlinearities. 

To  a  construct  a  Preisach  model  capable  of  characterizing  the  asymmetric  hysteresis  loop,  the  measure  /z(s) 
was  approximated  by  the  expansion  (4)  with  the  coefficients  ani  chosen  through  a  least  squares  fit  to  the  data. 
The  resulting  model  response  obtained  using  80  coefficients  accurately  characterizes  the  data  as  illustrated  in 
Figure  4.  Hence  the  Preisach  model  provides  a  technique  for  accurately  quantifying  the  nonlinear  constitutive 
properties  of  the  transducer  but  at  the  expense  of  a  large  number  of  nonphysical  coefficients. 


4.  Concluding  Remarks 

The  Preisach  techniques  employed  here  provide  a  general  methodology  for  characterizing  hysteresis  inherent 
to  the  piezoceramic  positioning  mechanisms  employed  in  atomic  force  microscopes.  While  illustrated  in  the 
context  of  a  stacked  actuator,  analogous  models  can  be  developed  for  a  cylindrical  transducer  by  incorporating 
the  Preisach  model  in  2-D  constitutive  relations  employed  in  general  thin  shell  theory  (see  [14]).  This  provides 
a  technique  for  constructing  a  variety  of  full  and  reduced-order  models  for  nanopositioners. 

The  advantages  provided  by  the  generality  of  the  technique  are  countered  by  the  large  number  of  nonphysical 
parameters  required  when  characterizing  biased  or  asymmetric  hysteresis  loops  along  with  the  extensions  to  the 
theory  required  to  accommodate  temperature  or  frequency  variations.  For  such  operating  regimes,  alternative 
models  of  the  type  described  in  [5,  10]  can  be  considered. 


x  10' 
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Figure  4.  Experimental  data  from  the  stacked  actuator  and  model  prediction. 
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